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ABSTRACT 



The spectrum of the normal Type la SN 1990N observed very early on (14 
days before B maximum) was analysed by Fisher et al. (1997), who showed that 
the large width and the unusual profile of the strong line near 6OOOA can be 
reproduced if the line is assumed to be due to C ll 6578, 6583A and if Carbon 
is located in a high velocity shell. This line is one of the characterising features 
of SNe la, and is usually thought to be due to Si ll. A Monte Carlo spectrum 
synthesis code was used to investigate this suggestion further. The result is 
that if a standard explosion model is used the mass enclosed in the shell at 
the required high velocity (25,000-35,000 km s^^) is too small to give rise to 
a strong C ll line. At the same time, removing Silicon has a negative effect 
on the synthetic spectrum at other wavelengths, and removing Carbon from 
the lower velocity regions near the photosphere makes it difficult to reproduce 
two weak lines which are naturally explained as C II, one of them being the 
line which Fisher et al. (1997) suggested is responsible for the strong 6OOOA 
feature. However, synthetic spectra confirm that although Si II can reproduce 
most of the observed 6OOOA line, the red wing of the line extends too far to 
be compatible with a Si ll origin, and that the fiat bottom of the line is also 
not easy to reproduce. The best fit is obtained for a normal SN la abundance 
mix at velocities near the photosphere (15,500-19,000 km s~^) and an outer 
Carbon-Silicon shell beyond 20,000 km s^^. This suggests that mixing is not 
complete in the outer ejecta of a SN la. Observations at even earlier epochs 
might reveal to what extent a Carbon shell is unmixed. 

Subject headings: supernovae: general - supernovae: SN 1990N - line: 
identification - line: formation - line: profiles 
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1. Introduction 

The distinguishing spectral feature of SNe la is the Si ll 6347, 6371A hne. This hne 
is strong in all photospheric-epoch spectra of SNe la, appearing near 6OOOA from before 
maximum to at least one month after maximum. Its strength suggested that a large 
fraction of the material ejected in a SN la explosion is in the form of Si and other so-called 
'Intermediate Mass Elements' (i.e. between CO and the Fe-group). These elements (mostly 
Si, S, and Ca) are synthesised in regions where thermonuclear burning of the progenitor 
Carbon-Oxygen white dwarf (WD) is incomplete, not leading to the production of ^^Ni 
(Nomoto et al. 1984). Since this is only possible if the burning wave that disrupts the WD 
is subsonic (aka a deflagration wave), the identification of this line was an important step 
in understanding that burning must must proceed subsonically in a significant fraction 
of the WD. Explosion models based on this assumption could successfully reproduce the 
light curve and spectra of SNe la starting from a Chandrasekhar mass CO WD (e.g. 
the defiagration model W7, Nomoto et al. 1984, Branch et al. 1985; or various delayed 
detonation and pulsational defiagration models presented by Woosley & Weaver (1994) and 
Woosley (1997). 

Fisher et al. (1997) (hereafter FBNB) re- analysed the earliest spectrum ever recorded 
of a SN la, that of SN 1990N on 26 June 1990, which is as early as 14 days before maximum 
B light (Leibundgut et al. 1991). Using a simple but powerful parametrised LTE code, 
FBNB reproduced the UV-optical spectrum satisfactorily, but noticed that the broad and 
fiat-bottomed profile of the 'Si ll' line, observed near 6040A, could not be reproduced using 
Si and an exponential radial dependence of the line optical depth. Other authors (Jeffery 
et al. 1992, Mazzali et al. 1993) encountered essentially the same problem, even though 
they used more sophisticated codes. As an alternative solution, FBNB showed that the line 
profile could be successfully modelled by a high velocity Carbon shell (26, 000 < v < 35, 000 
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km s^^). This proposed alternative identification is not in contradiction with the basic 
properties of successful explosion models, which usually have a layer of unburned progenitor 
material (CO in this case) at the top of the ejecta and hence expanding at the highest 
velocities. Such a layer would be most easily detectable for its effect on the spectrum when 
it is not yet very far removed from the momentary and receding photosphere, and therefore 
it is natural to expect that the very early spectrum of the normal SN la 1990N would be 
one of the best places to find evidence for it. 

This suggestion has very interesting consequences, both because it sets some limit to 
the mixing taking place in the highest velocity part of the ejecta and because it points out a 
possible source of error in measures of very early spectra of SNe la, especially if the spectra 
are of poorer quality than that of SN 1990N, a typical situation for SNe at high redshift. 

Independently of FBNB, we have been modelling the same UV-optical spectrum of 
SN 1990N in an effort to determine its epoch as accurately as possible so as to set some 
limits on the rise time to maximum, in view of the current debate about this parameter 
(Mazzali & Schmidt 2000, in prep.). We have used a Monte Carlo (MC) code based on 
that described in Mazzali & Lucy (1993), but modified to include an extended line list and 
photon branching (Lucy 1999, Mazzali 2000). Although we were trying to obtain a good 
overall fit of the spectrum, and did not pay much attention to small features, we noticed 
that in just about every synthetic spectrum we computed two rather weak features were 
always well reproduced as C II lines: a small absorption near 6350A, which is attributed to 
C II 6578, 6583A, and the absorption near 6900A, attributed to C ll 7231, 7236A with a 
contribution from O II 732lA to make the feature relatively broad. The strong absorption 
line near 6040A was matched by the Si ll line, but the line profile was not reproduced 
correctly, as in all previous work. In all models the abundances were homogeneous above 
the photosphere. The model we finally selected had an epoch t = 5 days (implying 
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t{Max) = 19 days), distance modulus /i = 32.00 mag, luminosity L = 2.30 ■ 10^^ erg and 
photospheric velocity Vph = 15750 km s^^. The synthetic spectrum is shown as the thin 
continuous line in Fig. 1. The velocity indicated by the two weak features in the red, if 
they are interpreted as C ll lines, is about 12000 km s~^, which is somewhat smaller than 
Vph- In fact the synthetic C II lines, especially that at 6350A, fall at a somewhat shorter 
wavelength than the observed ones. Nevertheless, the coincidence was striking. 

Thus we turned our attention to the FBNB paper. FBNB show fits to the entire 
spectrum using different compositions in their Fig.l, which was unfortunately printed too 
small to verify whether the two small features discussed above are reproduced as C ll. 
Nevertheless, in their Fig. 2, where they show two different ways to reproduce the 6040A 
feature - using Si and C, respectively, it is clear that neither the 6350A feature nor the one 
at 6900A are reproduced with the 'high velocity C shell' model. On the contrary, while in 
that model the C ll 6578, 6583A doublet gives rise to the strong and flat-bottomed 6040A 
absorption, a weaker, broad and also flat-bottomed absorption is also visible near 6500A, 
which is not present in the observed spectrum. This is exactly where the C ll 7321, 7236A 
doublet would fall if Carbon were located in a shell between 26,000 and 35,000 km s~^. 
The ratio of the equivalent widths of the two features is large, with a value of at least 
10, which is larger than the observed ratio but may result from the assumption of LTE 
since the redder doublet comes from a more highly excited level (16.33 v. 14.45 eV). We 
therefore considered it worthwhile to tackle once more the problem of what might give rise 
to the observed 6040A line using synthetic spectra. In the next sections we discuss various 
alternatives. 
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2. Can Carbon replace Silicon? 

We have shown that a fully mixed model reproduces two weak C ll features, and gives 
a less-than-satisfactory fit of the 6040A feature as Si ll. If we are to explain this feature 
as due to C II instead, the first necessary step is to remove Silicon from the mixture (as 
FBNB did in the central panel of their Fig.l) and ascertain what infiuence this has on the 
synthetic spectrum. In Fig.l we show two models obtained for a homogeneous composition: 
the dashed line is for a full W7 composition mix, including both Si and C, as discussed 
in the previous section. The peak at 3500A, which is not reproduced by our synthetic 
spectrum, is in the MMT data. Since the calibration at the edge of the spectrum is tricky, 
the poor fit at this wavelength does not necessarily mean that the model is wrong. 

The dotted line in Fig.l is a model where Silicon has been removed and its abundance 
(0.2 by mass) is assigned to Carbon instead. In this model both C ll lines, but especially 
the 6580A one, are too strong, and are not compatible with the observations. Removing 
Si, and not placing C in a high velocity shell, clearly destroys the fit of the 6040A feature. 
Furthermore, removing Si has negative consequences elsewhere in the spectrum, which 
cannot be remedied by introducing C at high velocities: 1) the weak Si ill 4553, 4568A 
line, observed near 4400A, a typical feature of all but the coolest SNe la before and near 
maximum, is now lost; 2) the strong absorption near 4800A has about equal contributions 
from Si II 5041, 5056A (falling near 4750A) and Fe ill 5156 + Fe ll 5169 A (falling near 
4850A) in the spectrum obtained for a W7 mix: when Si is removed, only the red part of the 
absorption is left, and this absorption is too weak to cause a strong re-emission peak, which 
is observed near 5100A; 3) the shape of the spectrum near 5500A is infiuenced mostly by 
Si III 5740A (near 5425A) and Si ll 5979A (near 5600A), although the broader absorption 
extending to 5200A is due to several lines of S ll. When Si is removed from the mixture, 
only the S II lines are left and the absorption in the model is too far to the blue. These 
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three points show clearly that Silicon cannot be completely removed from the relatively 
high velocity regions near the photosphere. But if Si is allowed in the mix, then the doublet 
6347, 6371A is the strongest Si line. Therefore, it appears unavoidable to conclude that 
Si II contributes significantly to the 6040A feature even at this early epoch. 

3. A Carbon-shell model? 

FBNB could fit the 6040A feature placing Carbon at high velocity, but since they 
used a parametrised model they did not quantify the mass of C required to obtain such a 
strong synthetic line. Since the C ll lines are weak even in a fully mixed model, this is an 
interesting numerical experiment. The model shown as a dotted line in Fig.l, where Si had 
been replaced by C, had a Carbon abundance of 0.25 by mass. Since the total mass above 
the photosphere in this model is 0.072Mq, the C mass was O.OISMq, and the strength of 
the C II 6578, 6583A line was comparable to that of the observed 6040A feature, but the 
velocity was wrong. The velocity can be reconciled if C is enhanced only for velocities 
larger than Vsh = 26000km s~^, as FBNB suggested. In our next synthetic spectrum we 
replaced Si with C only above this velocity, while between Vph and Vgh we just redistributed 
the Si abundance among all elements, proportionally to their respective abundances (the 
composition is dominated by Oxygen). Because of the steeply falling density {p ex r~^), 
when C is confined to very high velocities its mass is reduced. This reduces the C mass 
to only 1O~^M0. The result is shown in Fig. 2 (thin continuous line): both C II lines are 
so weak that they do not produce a noticeable feature in the spectrum. This calculation 
reproduces all essential features of the FBNB 'high velocity C shell' model, and at the same 
time is consistent with the hydrodynamical model of a SN la: the result is that not only 
the synthetic spectrum does not reproduce the strong 6040A feature, but the absence of 
Si creates the three problems listed in Sect. 2, and the displacement of C in velocity space 



means that the two weak features which are attributed to C ll are not reproduced either. 

Obviously, it must be possible to fit the 6040A feature as C ll by increasing the C ll 
line opacity in the high velocity shells. We show this as the dashed line in Fig. 2. The 
values by which we had to multiply the Sobolev optical depth to obtain this spectrum are 
an increasing function of velocity, ranging from 10 at 22500km s^^ to 10® at 30000km s~^, 
the outer limit of the assumed ejecta distribution. Although these values do give a very 
good fit to the observed feature, they really have no physical basis. If the increased optical 
depth is attributed to the number of absorbing C ll ions, this would lead to a completely 
unrealistic C mass, more than IMq, at these large velocities. Furthermore, all the problems 
raised by the absence of Si remain, and the two small features are not fitted at all. 

It is interesting to note that although we increased the optical depth for all C ll lines, 
we obtained a strong 6578, 6583A line, but not a strong 7231, 7236A line. This is because 
the latter doublet comes from a more highly excited level, whose population falls more 
steeply with radius, and so even larger factors would be required to enhance the line at very 
high velocities. 

This test confirms that the 6040A feature must be at least predominantly due to Si ll, 
and that placing Carbon in a high velocity shell does not give sufficient opacity, even if C 
dominates the composition. 

4. A Silicon model? 

We believe we have provided ample evidence that the 6040A feature cannot be entirely 
- or even mostly - due to C ll lines. As an alternative solution, can the distribution of 
Silicon be modified so that a better fit to the line can be obtained? If we look back 
at Fig.l we notice that the profile of the Si ll 6347, 6371 A line is too sharp, and that 
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absorption is missing both in the blue (between 5800 and 5900A) and in the red (between 
6100 and 6200A). These regions correspond to Si velocities of about 23000 and lOOOOkm/s, 
respectively. Therefore, we can expect that we can produce more blue absorption by 
increasing Si at high velocity, but the velocity required to fit the red part of the absorption 
is significantly smaller than that of the photosphere. 

In Fig. 3 we show a model where the Si abundance has been increased by a factor of 
3 above 23000km s^^, at the expense of all other elements. This gives a reasonably good 
fit to the blue side of the line. Other Si ll lines are weaker, so they are not much affected 
by this change, and the rest of the spectrum is not very different from that of Fig.l. As 
regards the red part of the line though, setting the photospheric velocity at 15750 km s~^, 
imposes the constraint that the reddest wavelength of Si ll is about 6050A. The velocity of 
the photosphere cannot be reduced to about 10000 km s^^ without changing the aspect of 
the synthetic spectrum completely (see Mazzali & Schmidt 2000, in prep.). Therefore we 
agree with FBNB that Si alone is unable to explain the observed width and profile of the 
6040A feature. 

5. A possible solution: Silicon and Carbon 

On the basis of the models presented above we suggest that Si is responsible for most 
of the 6040A feature in SN 1990N, and in particular for its extended blue side, but an 
alternative origin must be found for the red part of the line. At the same time, we have 
shown that the small feature near 6300A is most likely due to C ll at a near-photospheric 
velocity. Our line list, which was derived from that of Kurucz & Bell (1995) does not seem 
to offer any reasonably strong line with a rest wavelength of about 6450A, which would 
naturally explain the red extension of the 6040A feature by redshifting at the velocity of 
the photosphere. The nearest strong lines to the red of the Si ll doublet are indeed the 
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C II ones. A Carbon shell centred at about 20000km s^^ might explain the observations. 
The questions are can C be responsible for both the small absorption at 6300A (via 
its near-photospheric distribution) and the red wing at ~ 6100A (via an enhanced C 
abundance at v ~ 20000km s~^), and if so is the required C enhancement compatible with 
the hydrodinamical structure of a SN la? 

The answer is shown in Fig. 4. This is a model computed with the same parameters 
as the previous ones and the following distributions: Carbon: 0.05 by mass up to 

V = 18500km s^^, then increasing smoothly and reaching a peak value 0.50 above 

V = 21500km s~^; Silicon: 0.20 by mass up to t; = 21500km s~^, and then increasing to 0.5 
above v = 22500km s^^. The integration is limited to an outer velocity of 30000km s^^. 

The model reproduces the observed profile reasonably well. The observed flat bottom 
of the line may require more ad hoc adjustments of the distribution of the elements. 
The 'average' positions of the near-photospheric and of the high velocity components of 
the synthetic C ll 6578, 6583A absorption are marked in Fig.4. The near-photospheric 
component falls at a longer wavelength than the the observed 6300A feature. The difference 
in velocity is about 3600km s^^. If the adopted photospheric velocity is correct, this offest 
is a puzzle. We have no alternative identiflcation other than C ll to offer for that feature, 
since our line list offers only very weak Ti II lines with rest wavelength of about 6700A. 
But even if the 6300A feature is not C ll, our conclusion regarding the origin of the 6040A 
feature does not change. 

The scenario we are faced with is then the following: a mixed W7 composition (C=0.05 
and Si=0.20 by mass, respectively) holds out to v ~ 19000 kms, but outside that velocity a 
C-Si shell develops. The near-photospheric C abundance may be lower if the line at 6300A 
is not C II. With the exception of the Si shell, this is consistent with incomplete mixing in 
the SN ejecta. Maybe S ll lines are responsible for the extra absorption between 5800 and 
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5900A, but our line data do not support that. 

6. Conclusions 

We have shown that the strong and broad absorption near 6040A in the d -14 spectrum 
of SN 1990N must be predominantly due to Si ll 6347, 6371A. If we reproduce the line as 
C II, eliminating Si from the mixture as suggested by FBNB, other regions of the spectrum 
where Si lines are strong are not well reproduced. Also, the mass enclosed in the high 
velocity shell which could give rise to a blueshifted C ll line is very small, so that the 
synthetic line is much too weak, and the observed profile can only be reproduced if the 
optical depth of the C ll lines is increased by unrealistic factors (up to 10^). On the other 
hand, we have confirmed the result of FBNB that Si ll alone cannot reproduce the entire 
feature because the red edge of the line has a blueshift much smaller than that of the Si ll 
line at the velocity of the momentary photosphere. 

As a possible contribution for that part of the feature we also suggest C ll in a high 
velocity shell, but unlike FBNB we enhance C in a shell between 19000 and 30000km s~^, 
which contains a small (O.O4M0), but nevertheless large enough mass to give rise to a rather 
strong line without resorting to artificially large departure coefficients to increase the line 
optical depth. The increased C abundance in this outer shell does not significantly affect 
the total mass of C in the ejecta. At the same time, mixed C at near-photospheric velocities 
gives rise to two weak synthetic absorptions at 6250 and 6900A, which are a reasonable 
match for two observed features at 6300 and 6900A. The blending of the Si ll and C ll lines 
may give rise to the flat bottom of the observed 6040A feature, although our model does 
not reproduce that. 

The implications of our result are that an outer zone where Carbon is not fully mixed 
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does indeed appear to exist. This zone is narrow, and contains only a small mass, so 
the photosphere passes through it rapidly. Our models suggest that Si is also present 
there. Therefore, observing a SN la very early on, when the photosphere is still at 
V ~ 20000km s~^, would be very interesting. If both Si and C are present at high velocities, 
we do not expect that the feature should appear very different from what is observed in the 
spectrum of SN 1990N we have analysed in this paper. Alternatively, if only C were present, 
the 'Si ll' absorption would be much redder (A ~ 6150A) and caused mostly by C ll. This 
would require observations very soon after the explosion, roughly 16 days before maximum, 
which this may not be so very difficult with SNe at high redshift, whose observed evolution 
is slowed down by the factor {1 + z). Such observations would be extremely interesting 
because they would probe the outermost part of the ejecta and could thus clarify by direct 
observation to what extent mixing actually takes place. Finally, our findings must serve as 
a caveat against deriving SN la properties from the 'Si ll' line if the SN is observed very 
early on. 

It is a pleasure to thank Brian Schmidt for providing the observed spectrum in a form 
a modeller can easily work with. 
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Fig. 1. — The — 14d UV-optical spectrum of SN 1990N (thick sohd hne) is compared to: 
a) a synthetic spectrum computed with a W7-mixed composition (thin sohd hne), which 
reproduces the overaU spectrum, including several Si ll and Si ill features. The Si ll 6347, 
637lA line is too narrow, and two C II lines are seen: one, near 6250A, may be compatible 
with the feature at 6300A, while the other is a good match for the line at 6900A; b) a 
synthetic spectrum where Si has been replaced by C throughout the envelope (dotted line). 
The C II lines (especially the one at 6250A) are too strong, while many Si features are 
missing. 
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Fig. 2. — The — 14d UV-optical spectrum of SN 1990N (thick sohd hne) is compared to: a) a 
synthetic spectrum where Si is ehminated, but it is replaced by C only aX v > 26000km s~^, 
in analogy with FBNB (thin solid line). The C ll lines are too weak to leave a signature in the 
synthetic spectrum; b) a synthetic spectrum similar to the one above, but with artificially 
increased C ll line opacities. This fits the 6040A line, but fails to reproduce other Si ll and 
Si III lines and the two weak lines at 6300 and 6900A. 
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Fig. 3. — The — 14d UV-optical spectrum of SN 1990N (thick sohd hne) and a synthetic 
spectrum computed for an increased Si abundance above 23000km s"^(thin sohd hne). This 
reproduces the blue side of the 6040A hne, but not the red side. 
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Fig. 4. — The — 14d UV-optical spectrum of SN 1990N (thick sohd hne) and a synthetic 
spectrum computed with increased Si and C abundances above about 20000km s~^ (see text 
for details). This reproduces most observed features, including the broad 6040A line and 
the two weak C ll lines, although the near-photospheric component of the 6578, 6583A line 
does not match the wavelength of the observed 6300A absorption. 



